Missing physics in stick-slip dynamics of a model for peeling of an adhesive tape 
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It is now known that the equations of motion for the contact point during peeling of an adhesive 
tape mounted on a roll introduced earlier are singular and do not support dynamical jumps across 
the two stable branches of the peel force function. By including the kinetic energy of the tape in the 
Lagrangian, we derive equations of motion that support stick-slip jumps as a natural consequence 
of the inherent dynamics. In the low mass limit, these equations reproduce solutions obtained using 
a differential-algebraic algorithm introduced for the earlier equations. Our analysis also shows that 
mass of the ribbon has a strong influence on the nature of the dynamics. 

PACS numbers: 05.45.Pq, 62.20.Mk, 68.35.Np 
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It is well established that there are three different 
modes of failure during peeling of an adhesive tape from a 
substrate 0. At low applied velocities, the peeling front 
keeps pace with the pull velocity and the failure mode 
is cohesive. At high pull velocities, the failure is adhe- 
sive. In between these two regimes, there is an intermit- 
tent mode of failure corresponding to stick-slip dynamics 
0,0 accompanied by a characteristic audible noise. The 
stick-slip nature suggests that this regime is unstable. 
Indeed, the strain energy release rate exhibits two sta- 
ble branches separated by an unstable branch. Detailed 
studies by Maugis and Barquins Q and others 0, 0, 
show that the nature of the wave forms of the pull force 
exhibits sinusoidal, sawtooth and highly irregular wave 
patterns. Gandur et al. |(| have carried out a dynamical 
analysis of the force waveforms and of acoustic emission 
signals, and report chaotic dynamics at the upper end 
of the pull velocities . Analysis of acoustic signals has 
also been carried out by Ciccotti et al. [U- 

Stick-slip behavior is commonly observed in a num- 
ber of driven systems such as the Portevin-Le Chatelier 
(PLC) effect @, frictional sliding 8], and earthquake 
dynamics 0. Stick-slip is characterized by the system 
spending most part of the time in the stuck state and 
a short time in the slip state. One common feature of 
such systems is that the force exhibits " negative flow rate 
characteristic" (NFRC). In fact, studies of dynamics of 
such systems, including that of the adhesive tape, use 
the macroscopic phenomenological NFRC feature as an 
input, although the unstable region is not accessible. 

Maugis and Barquins were the first to look at the 
peeling problem from a dynamical angle. Later, a de- 
tailed study was carried out by Hong and Yue [3( on a 
three variable model originally introduced by Maugis [lj} 
using an "N" shaped function that mimics the peel force 
function. They report that the system of equations dis- 
played periodic and chaotic stick-slip solutions. However, 
it was later recognized that the stick-slip oscillations were 
not obtained as a natural consequence of the equations 
of motion 0, 0, 0] as the jumps in the rupture speed 
were introduced externally once the velocity exceeded the 
limit of stability. Later, Ciccotti et al. pj| interpret 
stick-slip jumps as catastrophes with the belief that the 



jumps in the rupture velocity cannot be obtained from 
these equations of motion [J] . Recently we derived these 
equations starting from a Lagrangian and showed |Tl) 
that these equations are singular and fall in the category 
of differential-algebraic equations (DAE) (Lj| requiring a 
special algorithm. Using a DAE algorithm, we showed 
that stick-slip jumps across the two branches arise in a 
pure dynamical way. The dynamics was also shown to 
be much richer than anticipated earlier. 

However, even as the DAE algorithm offers a mathe- 
matical framework for obtaining solutions for these sin- 
gular equations, it is difficult to provide any physical 
interpretation for the "mass matrix" that removes the 
singularity. Thus, a proper identification of the missing 
physics responsible for the absence of dynamical jumps in 
these equations still remains to be addressed. As we shall 
comment later, this is a necessary step for understand- 
ing the origin of acoustic emission (AE) during peeling, a 
problem that has remained unresolved. A reexamination 
of the earlier derivation showed that the kinetic energy 
of the stretched part of the tape was ignored. Here, we 
show that the inclusion of this additional kinetic energy 
removes the singularity thus converting them into set of 
ordinary differential equations (ODE). Further, stick-slip 
jumps emerge as a natural consequence of the inherent 
dynamics itself. Apart from reproducing the DAE solu- 
tions for low mass limit, our analysis shows that the mass 
of the tape has a considerable influence on the nature of 
the dynamics. 

A schematic of the experimental setup is shown in Fig. 
ID^a). An adhesive roller tape of radius R is mounted on 
an axis passing through O normal to the paper and is 
driven at a constant speed V by a couple meter motor 
positioned at O' . The pull force F acting along line join- 
ing the contact point P and O' subtends an angle 9 to 
the tangent at the point P. As the contact point P mov- 
ing with a local velocity v can undergo rapid bursts of 
velocity during rupture, the peeled length of the ribbon, 
L = PO' is not fixed. We note here that the peel force 
often called the force of adhesion denoted by / is what is 
measured in experiments in steady state conditions. Let 
the distance from the center of the roller tape O to the 
motor O' be I and a the angle subtended by PO with 
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FIG. 1: (a) Schematic plot of experimental setup, (b) Plots 
of f(v, V) for V = 1 (solid curve), V — 4 (dotted curve). 
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FIG. 2: (a), (b) The v — F phase plots for I = W~ 5 ,V = 1 
corresponding to the DAE and ODE solutions for m = 10 -4 . 
Solid line shows f(v, 1). 



the horizontal line OO' . Let I be the moment of inertia 
of the roller tape, w the angular velocity, k the elastic 
constant of the tape, and u the elastic displacement of 
the tape. As the contact point is not fixed, u> = a + v/R. 
The geometry of the setup gives L cos 9 — —I sin a and 
L s'm9 = I cos a — R. As the peeling point P moves, the 
pull velocity V is the sum of three contributions Q, i.e., 
V = v + u — L, which gives 



v = V — ii + L = V — u — R cos 9 a. 



(1) 



Then, the Lagrangian is C = Uk — Up with the kinetic 



energy given by Uk = klw 



and the potential 



energy by Up — hku 2 . [\lu 2 is the kinetic energy of the 
roller tape and |mi 2 ignored earlier, arises due to the 
kinetic energy of the stretched part of the tape, u refers 
to the time derivative of u. In principle m (m = moL(t) 
where m is the mass per unit length) depends on time 
through L. However, for all practical purposes m can be 
treated as constant.] We write the dissipation function 
as 1Z = <5(u, V) = J f(v, V)dv, where f(v, V) assumed to 
be derivable from a potential function V), physically 
represents the peel force, and is taken to depend on the 
pull velocity and rupture speed as in 

Using the Lagrange equations of motion, ^ (§§) — 
|£ + P = and 4 (SS) - 2£ + |5 = 0, and using 
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(a, d, u, u) as generalized coordinate, we get 



v R cos ( 



:f(v,V), 



(2) 



R I (1- cos 

' -[f(v,V)-ku(l-cos9)}. (3) 
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Note that the right hand side of Eq. (PJ) is the algebraic 
constraint in Eq. (10) of Ref. [13 or Eq. 5(d) of Ref. 
0. These equations in their present form are still not 
suitable for further analysis as they have to satisfy the 
constraint Eq. (JJJ. In the spirit of classical mechanics 
of systems with constraints (see Ref. 53 )j we derive the 
equation for the acceleration variable v in the constraint 
equation by differentiating Eq. and using Eqs. (J5J), 
and ©. This gives equations of motion for a, u, u and v 
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(4) 
(5) 



u = V — v — R cos 6 a, 

v = — u + R sin 9 9 a — R cos 8 a. 



(6) 
(7) 



Using Eqs. ©, © in Eq. Q, we get 
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We retain the dynamization scheme introduced earlier 
[Tl| adopted from the PLC effect [l5j wherein the differ- 
ence between the maximum and minimum of f(v, V) is 
assumed to decrease with increasing V. The parameter- 
ized form of f(v, V) used here is given by 



f(v,V) = 402t- - 34 + 171 V °- 16 + eSe^ 7 - 7 ) - 2U 1 - 5 

-(415-45U a4 -0.35U 215 )w - 5 . (9) 



Equation @ mimics the general trend of the experimen- 
tal peel force function and is essentially the same as used 
earlier jlll ] except that it accommodates larger excur- 
sions of trajectories arising from the introduction of the 
additional time scale [Fig. IHb)]. The fixed point of Eqs. 

El El El given by a = 0,u = V/R,u = f(V,V)/k 
and v — V becomes unstable when V is such that 
/'(V, V) < leading to a Hopf bifurcation. 

We have solved Eqs. (@J, JSJ, and JSJ) by adap- 
tive step size stiff differential equations solver (MATLAB 
package) and studied the dynamics over a wide range of 
values of / (kg m 2 ), V (m s _1 ) for m (kg) ranging from 
1CT 4 to 0.1 (keeping k = 1000 N m" 1 , R = 0.1 m, and 

1 = 1 m). (Henceforth, we suppress the units for the 
sake of brevity.) Here, we present results (obtained af- 
ter discarding the initial transients) for a few represen- 
tative values of the m when (/, V) are at low and high 
values. We show that for low m, we essentially recover 
the DAE solutions reported earlier Further we show 
that mass of the tape has a strong influence on the nature 
of the dynamics. 

A rough idea of the nature of the dynamics can be 
obtained by comparing the frequency Q u = (fc/m) 1 / 2 as- 
sociated with u with Q a — (Rf jl) 1 ! 2 corresponding to 
a. Since / (in N) is limited to 180-280, the range of 
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FIG. 3: (a), (b) Phase plots of v - F obtained for m = 10" 4 
and 0.1 for I = 10~ 2 , V = 1. Solid line shows f(v, 1). (c), (d) 
Plots of v(t) for m = 10~ 4 and 0.1. (e) The largest Lyapunov 
exponent, A for m — 0.1, / = 10~ 2 , V = 1. (Units of v, V are 
in m s _1 , F in N, / in kg m 2 , t in s and A in s -1 .) 



n a (s^ 1 ) is 1342-1673 for small I (~ 10" 5 ) decreasing 
to 42-53 for large I (0.01). In comparison, fl u (s _1 ) is 
3162 for rn = 10~ 4 decreasing to 100 for m = 0.1. Thus, 
one expects the nature of solutions to be influenced with 
increasing m for fixed I which also depends of V ^l| . 

Here, we present our numerical results. Consider the 
low mass limit, i.e., low kinetic energy of the tape. Then 
the RHS of Eq. J2J is small, and turns f(v, V) — ku(l — 
cos 9) f» f(v, V) — F (1 — cos 6) = (which is the algebraic 
constraint that makes the equations singular). Thus, one 
expects that the DAE solutions are reproduced for small 
m which we have verified for the entire range of values of 
I and V studied previously (For numerical calcula- 
tions, we have used m — 10~ 4 as the low mass limit.) As 
an example Figs. Ufa) and[2Ib) show the phase plots in 
the v — F plane obtained using the DAE algorithm and 
ODE Eqs. Q, ©, ©, and © for m = 10~ 4 respec- 
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FIG. 4: (a), (b) Plot of v(t) for m = 10~ 4 and 0.1 for 
/ = W~ 2 ,V = 4. (c), (d) The corresponding phase space 
trajectories. f(v, V) is shown by a bold line. (Units of v, V 
are in m s~\ F in N, / in kg m 2 and t in s.) 



tively keeping 7=10 5 ,V = 1. It is evident that the 
DAE solution is similar to the ODE solution. 

Much more complex dynamics emerges as a result of a 
competition between this additional time scale and other 
time scales present in the system. Consider the results for 
m = 10~ 4 and 0.1, for I = 10~ 2 and V = 1. The small 
m plots are provided for the sake of comparison as they 
essentially correspond to the DAE solution. Consider 
the phase plots in v — F plane shown in Figs. OJa) for 
m = 10~ 4 and Fig. IHJb) for m = 0.1. It is clear that the 
influence of increasing m is considerable. In particular, 
note that the sharp changes in F in the v — F plot [Fig. 
Efa)] at the upper end of v for small m are rendered 
smooth for large m case [Fig. E^b)]- Indeed, the effect of 
the additional time scale due to finite mass of the tape 
is also evident in the plots of v(t) for the low and high 
mass cases shown in Figs. |3[c) and Eld) respectively. 
Finally, it is clear that the phase plot [Fig. Efb)] fills the 
space and is suggestive of chaotic dynamics. The chaotic 
nature can be ascertained by calculating the Lyapunov 
spectrum. Using the QR decomposition method [l6| , we 
have calculated the Lyapunov spectrum and find a large 
positive exponent with a value ~ 4.4s" 1 [Fig. EJe)]. 

Increasing m does not always increase the level of com- 
plexity of the solutions. As an example, Figs. 0a) and 
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FIG. 5: Phase diagram in the I — m plane for various values 
of V. Periodic •, long periodic x, chaotic (marginal) o and 
chaotic *. (Unit of / is in kg m 2 and m in kg.) 

01b) show plots of v(t) for m = 10~ 4 and 0.1 respectively 
for / = 10~ 2 for V = 4. While the solution for small m is 
similar to that of DAE which exhibits several sharp spikes 
in velocity [Fig. El a )]> for large mass (m = 0.1), v(t) is 
surprisingly simple and is periodic [Fig. EJb)]. Indeed, 
this is better seen in the phase plots v — F for m = 10~ 4 
and 0.1 shown in Figs. 0{c) and ^d) respectively. In 
contrast to the low mass v — F plot which is chaotic (see 
similar DAE solution in Fig. 5(c) in Ref. llj), a simple 
limit cycle emerges for m = 0.1. As the nature of the 
dynamics can vary from a simple limit cycle to a chaotic 
attractor as the three parameters are varied, these re- 
sults can be summarized as phase-diagrams in the / — m 
plane for different values of pull velocities, V, as shown 
in Fig. JSJ. Apart from the chaotic state (*) seen for a 
few values of the parameters, for most values, the system 
is periodic (•) and a few other values, the attractor is 
long periodic (x). We also find a (marginally) chaotic 
attractor (o) for V = l,m = 10~ 4 ,/ = 10~ 3 for which 
the positive Lyapunov exponent is ~ 0.03 s" 1 (which is 
much beyond the error in computation). 

In summary, we have demonstrated that the missing 
time scale arising from the kinetic energy of the stretched 
part of the tape plays an important role in the peeling 
dynamics of the adhesive tape. As the inclusion of this 
term lifts the singularity in the equations of motion hith- 



erto considered, stick-slip jumps across the two resistive 
branches emerge as a consequence of the inherent dy- 
namics. Further, our study shows that the mass of the 
tape has a strong influence on the nature of the dynam- 
ics. For low pull velocities, and high /, the complexity 
increases, i.e., trajectories that are not chaotic for low 
mass become chaotic with increasing m. In contrast, for 
high V, the trajectories that are chaotic for low m are 
rendered nonchaotic with increase of m. 

Apart from resolving the central issue, the inclusion of 
kinetic energy of the tape provides a mechanism for con- 
verting the potential energy stored in the stretched tape 
into kinetic energy and hence provides a basis for explain- 
ing acoustic emission during peeling. This involves first 
extending the model to include the spatial degrees of free- 
dom of the peel front and including an additional (rup- 
ture) velocity dependent dissipation function to mimic 
the AE energy dissipated (along the lines in 17]). The 
extended model also helps to analyse the contact line dy- 
namics of the peeling front which in itself has not been 
addressed so far. Preliminary results 01 show that the 
energy dissipated occurs in bursts similar to the nature of 
AE signals seen in experiments p| . Here, it is worthwhile 
to comment on the dynamization of the friction law. The 
physical origin of this can be attributed to the viscoelas- 
tic nature of the fluid which in turn implies frequency 
dependent elastic constant. Thus as higher pull speed 
allows lesser time for internal relaxation to be complete, 
the viscoelastic fluid behaves much like an elastic solid. 
Clearly, a rigorous derivation of the peel force function 
from microscopic considerations that includes the effect 
of the viscoelastic glue at the contact point is needed to 
understand the dynamics appropriately. 
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